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Landscapes that are rhythmically dissected by natural drainage channels exist in various 
geologic and climatic settings. Such landscapes are characterized by a length-scale for 
the lateral spacing between channels. We observe a small-scale version of this process 
in the form of beach rills and reproduce channelization in a table-top seepage experi- 
ment. On the beach as well as in the experiments, channels are spontaneously incised 
by surface flow, but once initiated, they grow due to water emerging from underground. 
Field observation and experiment suggest the process can be described in terms of flow 
through a homogeneous porous medium with a freely shaped water table. According to 
this theory, small deformations of the underground water table amplify the flux into 
the channel and lead to further growth, a phenomenon we call "Wentworth instability" . 
Piracy of groundwater can occur over distances much larger than the channel width. 
Channel spacing coarsens with time, until channels reach their maximum length. 



1. Introduction 

Nature frequently exhibits regularly spaced channels that are formed by fluvial erosion. 
Figure n shows one such example with approximately regularly spaced drainage outlets. 
The spacing in naturally occurring examples varies from tens of centimeters to tens of 
kilometers, spanning five orders of magnitude. Some examples are more nearly periodic 
than others, but in either case there is a typical length-scale for the lateral separation 
of streams. This length-scale is of interest not only because it determines the shape 
of landforms around us, but also because it se ts one of t he scales of variation in the 
heterogeneities of sedimentary deposits ( Tailing et a^Jll997| . 

The ubiquity of this phenomenon further motivates its study. Periodically spaced 
drainage outlets are observed in various climates and geologic settings, such as mou ntain 
belts (Hovius 199^ faul t blocks ([Tailing et al. 1997), submarine canyons ( Inman e^l\ 



19791 lOrange era/]ll994l). valleys of volcanic islands (IWentworthlll928t iKochel fc Piped 
1986*), sapping canvons (jHoward et a/.lll988l: ISchumm et aJll995|) . and Martian gullies 



|Malin & Edgett 2000). The ubiquity of the periodicity implies a common phenomenon 
but not a common cause. Indeed, as we now proceed to review, theoretical formulations 
of erosion initiated by seepage and overland flow indicate that characteristic spacings 
are expected, but not necessarily for the same reasons. It is unclear whether the various 
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Figure 1. Elevation map of sapping canyons, Okaloosa, Florida jSchumm et o/. 1995"). This is 
one example of periodically spaced channels in nature. The image is 36km wide and the spacing 
is roughly 5km. 



forms of periodic channel spacing are due to a few basically similar mechanisms or if there 
are many essentially different explanations. There are many conceivable explanations for 
the periodicity. We approach this rich phenomenology by a detailed study of one specific 
example that is particularly accessible to study. 

Broadly stated, two "end-members" delimit processes of channelization llKirkbv fc ChorlevI 
ll967tlDunnjll98(]l) . At one extreme is erosion due to overland flow, wherein shear stresses 
imposed by a thin sheet of water flowing downhill act to erode a surface ('Horton"l945'). 
At the othe r extreme is erosion due to groundwater flow ( Kirkbv fc Chorlev 1967; Dunn^ 
1198(11199(1) . When a subsurface flow erodes a porous material as it emerges from the 
material, the process is called seepage erosion (or sapping). Overland flow and seepage 
erosion are illustrated schematically in Figure [3 

A co ntinuum model for overland flow was given in a landmark paper bv lSmith fc Brethertonl 
It was shown that under certain circumstances the eroding surface is unstable to 
sinusoidal channel-like patterns. However they found that the most unstable wavenumber 
was infinite, implying that t he most unstable wave length is at the microscopic granular 
scale. Noting this deficiencv. 'i zumi fc Parker! l) 1995 1) provided an alternative theory based 
on two components: a threshold condition for erosion and analysis of the full equations 
for shallow overland flow. The inclusion of the threshold condition yields an instability 
at a finite wavelength, a finite distance from the divide. This result shows that periodic 
channelization can occur entirely as a consequence of its own dynamics. Under differ- 
ent circumstances, a fi nite wavelength can even emerge without any threshold condition 
fizumi fc Parkedl200(]() . 

The theory of seepage erosion is less well developed. As stated by iDunnd l)l99(]|) . 
"there have been no formal quantitative studies of the relation between the hydrogeologic 
propert ies [of material undergoing seepage erosion] and the average spacing of channels." 
iDunnd (jl980, 1990), however, provides a physically appealing description. The pressure 
head in a horizontal or tilted plane leads to high pressure gradients at places of high 
local curvature, and hence presumably to enhanced erosion. This instability leads to 
the growth of channels. A similar but more complex situation is depicted in Figure 
where a channel on an inclined slope is eroded below the point where the water table 
intersects the surface. Once a channel forms, internal flow lines are focused toward it. 
As we show below, the focusing occurs over a large but finite transverse extent, which 
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Figure 2. Schematic view of two erosive processes: overland flow (a) and seepage erosion (b). 
In overland flow a thin sheet of water (gray) flows over the surface, eroding it as a result of the 
shear stresses induced on the surface by the flow. In seepage erosion the fluid flows through the 
porous subsurface and erodes surface material after it emerges at the surface. Both processes 
can produce channels (the thick black lines). 



is presumably related to a characteristic channel spacing, the presence of which requires 

no heterogeneities in material properties. 

Documented experimental studies go back at least a century JjaggaiJIlQOSi). In the 
1980 ' s, se epage exper iments llKochel eit al\ Il985l: iKochel fc PipeilllQSfit ISchumm et a.l\ 
iHow ard 1988; Kochel et QiJll988il were frequently concerned with analogy to Mar- 
tian valleys. The question of whether ancient martia n channels are created by seepage 
or vast amounts of surface water re mains open today llBakerlll982l iLaitv fc Malinlll98iTl: 
lBakeilll99nt lAharo nson et a^Jl2002l) . Our own experiments are motivated by terrestrial 
analogs. Previous experimental work suggests that groundwate r piracy plays a dominant 
role in the development of channels. iHoward fc Mc 1/anel lll 988l) studied the rate of seep- 
age erosion in a narrow two-dimensional flow tank. iHowar dl (|l994lT99l has carried out 



computer simulation modeling of valley development by groundwater sapping. 

The present study is a first report of our effort to conduct field studies, experiments, 
and theory simultaneously to address the formation of channels, and rhythmic channel- 
ization in particular, by fluvial erosion when substantial groundwater flow is present. The 
following three sections cover, respectively, field study, experiment, and theory. The last 
section contains conclusions that we are able to draw so far. 



2. Field Study: Beach Rills 

Beach rills can ser ve as inspirati on for larger-scale drainage channels 

forming over geologic time periods (|HigginFJIl98l . While incised channels are not un- 
common on sand beaches, periodic channelization is rarer. Two suitable sites where small 
incised channels form when the tide recedes are located on Cape Cod. One is at Province- 
town, Massachusetts, a few minutes walk north of the main pier, the other at Mayo Beach 
in Wellfleet, Massachusetts. At high tide the sand is soaked with water, which seeps out 
as the water recedes. The tidal difference at both locations is typically 3 meters. All chan- 
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Figure 3. Beach rills observed in Provincetown, Massachusetts on 30 August 2001. The 
shoreline is parallel to and below the picture. The distance from left-to-right is approximately 

2 m. 

nels are erased when the tide rises again and they reform at different locations during 
the next tidal retreat. One can observe the generation of the channels on a daily basis. 

The channels found on the beach are sharply incised (Figure They have steep side 
walls so that their cross section becomes almost rectangular and they are about one cen- 
timeter wide. The sand grains are fairly distributed in size, around a mean of about 1mm. 
At the Provincetown site, the average channel spacing, based on over 200 measurements, 
is about 40cm. The spacing is measured near the head of the rills. Depending on the 
height of the tide, we observe 1-5 generations (bands) of channels, with typical lengths of 
0.7-2m. Many channels begin at pebbles, others start without visible perturbation. Fur- 
ther downhill, channels often merge and braid. The beach slope is 10-12%, terminated 
by a 1% slope that emerges only at low tide. At the Wellfleet site, the spacing is about 
Im, based on 38 channels, and the slope is 4%, also terminated by a lower slope. There 
are no pebbles at this site. 

We have observed channels several times during their initiation and growth at the 
Provincetown site. A vertical test hole allows us to monitor the ground water level as 
a function of time. The water table inside the beach sinks 1cm every 3 minutes dur- 
ing the hours before channel initiation. During the same time, the sea level sinks about 
three times as fast. Hence, the height difference between water table and sea level in- 
creases steadily until, about three hours after high tide, the pressure head suffices to form 
channels. 

The channels clearly form outside the water and are not directly related to the waves 
that hit the beach. Waves arrive occasionally and wash water up the beach, subsequently 
leading to a sheet of water flowing downhill. Pebbles interfere with this flow and tem- 
porarily cause cone-like flow patterns. These disturbances in the flow can lead to an 
incised streak downhill of the pebble which eventually grows into a channel. The initial 
incision happens within seconds or less and leaves a streak on the order of 10cm long. 
Not all channels are initiated this way. Another frequent initiation mechanism results 
from "fingers" of runoff water from any sources higher up (water ponding behind stones 
or from higher- lying, previous channels). Areas without surface defects sometimes do not 
develop any channels. At other locations even tiny perturbations suffice for initiation 
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Figure 4. An incised channel on the beach filled with running water. The image is about 
15cm wide. A pebble is located at the channel head. The water in the channel emerges from 
underground. 

and empty areas are "filled in" with channels. As far as one can see on the surface, the 
periodicity is established promptly (within a minute or so). 

Once a channel has initiated it grows downhill. Sand is washed away at the lower end of 
the channel by water running in the channel. At the downward end of the channel, where 
the longitudinal extension occurs, the water appears to be fully loaded with sediment. 
There is mass wasting on the sides so that the channels widen slowly; they also deepen 
slowly. Sand grains are transported by the water flowing in the channel. The rills grow 
to their final size within minutes to tens of minutes. 

Figure 01 shows a channel filled with water. The channel could be fed either by surface 
water that seeped out above the channel, or directly from groundwater provided from the 
bottom and sides of the channel. We found that they are directly fed by groundwater, for 
the following reasons. The film of water above the channel shows no movement visible to 
the eye. Small grains, placed by us on the surface, diffuse only slowly at speeds that are 
insufficient to provide the large amounts of water seen inside the channels. The channels 
carry water for quite some time after the initiation and long channels continue to be filled 
with water when the surface at its head is already dry. All this shows that the channels 
are fed directly by groundwater. 

The beach rills provide a natural example of rhythmic channelization and the obser- 
vations suggest that they require seepage for their explanation. 



3. Seepage Experiments 

3.1. Experimental Setup and Channels 

Prior experimental studies of seep age erosion llKochel et oiJIlOS^HKochel fc PipeilllQSfil: 
iHowardlllQSSt iKochel et aZ.lll988 ll have observed a few periodic channels. Our own ex- 
perimental setup is shown in Figure [S] The plexiglass box is 120cm wide and effectively 
90cm long. It is filled with 0.5mm glass beads. Water is pooled in the back and kept at a 
constant level using an overflow tube. The water seeps through the granular medium and 
creates channels on the slope, observed by a CCD camera. Sideways illumination makes 
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Figure 5. The experimental setup, a) A two-dimensional schematic cross section, b) The 

actual apparatus. 



the channels clearly visible in the image. The experiments are carried out with a steady 
water reservoir, unlike the beach where the water table is receding. 

A ledge protrudes at the bottom to increase the flux through the system. A minimum 
flux is required to dislodge a grain on an inclined surface which is below its angle of 
repose. The ledge makes it possible to achieve the threshold in a finite sized experimental 
apparatus, by allowing us to increase the height of the water column in the reservoir above 
the height of the granular bed. The resulting flux is also similar to that observed at the 
beach when channels are formed. Thus the ledge allows us to reproduce the processes 
observed on the beach in a finite sized experimental system. We did not observe any 
channels without the ledge. 

The capillary height in the medium is measured to be 25mm. The permeability is 
determined from flow through a U-tube. A water column of height A/i seeping through 
a porous medium over a distance L flows with a velocity of v = k/S.h/ L. The seepage 
coefficient is 3 mm/s. Sand from the beach, with its moderately larger grain size, has 
a moderately higher permeability, k ~8 mm/s. 

Before the beginning of each experiment a planar, sloped surface is created, then the 
water column in the back is filled to a prescribed height hw^ and the camera is turned on. 
The height hw of the water column is measured from the bottom of the tank. Figure 
shows examples of laboratory channels. These images are taken less than half an hour 
into the experiment. FigureEt shows about 10 drainage channels. FigureEli corresponds 
to a slump observed at high slopes. In this paper, we only report on patterns of fairly 




Figure 6. The outcome of four laboratory experiments done at various slopes s and water 
column heights hw- The water source is toward the top side of each image. Illumination is 
from both sides. The image width is 120cm. (a-c) Examples of multiple drainage channels, (d) 
Example of a slump observed at high slope. 



straight channels. Other erosion patterns, some of which are found bv lDaerr et all (|2003() 
in a different experimental setup, are observed in various parameter regions. 

Channels such as those in FigureEtL-c are observed within a limited range of slopes and 
water column heights. Figure summarizes the parameter regions for which experiments 
were carried out. In our early experiments, parameters were not accurately controlled, but 
Figure[7|corresponds to a later set of experiments conducted with a constant experimental 
setup. Channel formation is observed for slopes <17%. This slope is much less than the 
angle of repose of completely dry or saturated glass beads which is approximately 43%. 
For the particular level of box filling, channels are not observed below a water column 
height of 12cm. 

The mere observation of regularly spaced channels in the experiment immediately 
constrains possible explanations. The medium is composed of mono-disperse, cohesionless 
glass beads. The fluid is clean water. This is all that is required to create periodic channels. 
Reproducing channelization in a clean and controlled environment is a major step in 
studying the phenomenon and provides the basis for theoretical modeling. 

3.2. Temporal Evolution 

Figure IHl shows a time sequence of six images of the same experimental run. The channels 
are initiated after a sheet of water flows down the surface. The sheet disappears and the 
channels deepen, widen, and grow lengthward. This is similar to the beach rills, where 
channels are also initiated by surface water, then grow downhill by washing away sand 
at the lower end of the channel. Unlike the beach, the experimental channels, in addition 
to growing downward, also grow headwards until they stop at the point where the ledge 
ends. Channels grow initially (quickly) downwards and then slowly upwards (headward). 
After the channels have reached their maximum length they still widen and occasionally 
meander and fan at the lower end. 



N. Schorghofer, B. Jensen, A. Kudrolli, and D.H. Rothman 





18r 




17 








16 


-1—' 






15 


■(1) 


I 




C 

E 


14 






o 


13 


O 


i_ 




0) 


12 


■5 






11 













Channels 



••• 



No Channels 



10 12 14 



• 


Slumping 


• • 




□ □ 




□ 




□ □ 


16 


18 20 



Grade (%) 



22 24 



Figure 7. Observed behavior as a function of slope s and water column height hw- Filled circles 
(•) correspond to spontaneous channelization. High slopes lead to slumping of the material (□). 
No channels (0) form at low pressures and slopes. 




Figure 8. Periodic channels observed in the experiment. The water source is toward the top 
side of each image, a) An initial sheet of water, b) still no visible channels, c) incision of channels 
(marked with arrows), d) channels grow, some additional ones have formed, and others have 
dried up, e) continued growth, and f) the channels after about half an hour. New channels, not 
visible in the previous image, are marked with arrows. Channels that have discontinued growth 
are marked with crosses. 
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Figure 9. The position of major channels in two experimental runs as a function of time. Only 
some of the initial channels grow into large channels and coarsening occurs with time. The final 
spacing is more periodic than the initial state. Part (a) corresponds to the experiment shown in 
Figure H 



The spacing of initial channels, Figure|S|;, is al so apparently n ot random. It may result 
from an instability in the initial sheet flow (e.g. lHuDPerjll982() or from the interaction 
of the sheet flow with the eroding surface. Any threshold for incision is likely related to 
the sheet flow, rather than to the force exerted on the grains by water emerging directly 
from underground. 

Figures IHI and ini show spatial coarsening oyer the period of channel deyelopment. Only 
some of the initial channels grow into larger channels, while others cease their growth. 
For example, in Figure |Ht there are seyeral smaller channels yisible that haye stopped 
growing. This is especially the case when another channel is nearby. Growing channels 
tap into the groundwater supply of their neighbors, which then cease to grow due to the 
lack of groundwater. On the other hand, new channels still form after the initial incisions, 
especially when far from existing channels. Figure shows two examples. Both of these 
processes - the competition of nearby channels and the formation of new channels away 
from existing ones - eventually lead to a more regular spacing. 

In experiment (a) of Figure El there are 11 initial channels, of which 5 eventually cease 
to grow. The other 6 initial channels grow into major flnal channels along with 2 channels 
that formed at a later stage. In experiment (b) there are 11 initial channels, of which 3 
cease to grow and two additional channels form at a later time. These numbers indicate 
that many of the big flnal channels have grown directly from small initial channels. How- 
ever, the selection of which channels dry up and where new channels form is influenced 
by the distance to neighboring channels, so that the spacing is dynamically created or 
at least dynamically grown from the initial state. Channel locations have more of an 
imprint of the initial state than does the spacing between them. Incision behavior and 
groundwater piracy contribute to the final state. With more separation in time and space 
between initial and flnal conflguration one would expect an even larger spacing, which 
has little to do with the initial state. 



4. Theory and Numerical Simulation 

4.1. Governing Equation 

As a result of the observations on the beach and in the experiment, the problem can 
be formulated theoretically in terms of flow through a porous medium. After channel 
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Figure 10. Boundary conditions for a slope on top of an impermeable layer. The velocity 
component normal to the surface is denoted by v^. For the beach, the bottom boundary must 
be replaced by tailwater. The figure also defines ho, the height of the intersection between water 
table and surface, and vt, the absolute value of the velocity at this intersection. 

initiation, the supply of water to the channel is governed by groundwater flow. The 
erosion at the channel boundaries, especially at the lower end of the channel, is governed 
by the amount of groundwater collected in the channel. Our aim is to determine the 
extent of the groundwater piracy which limits the growth of channels. 

The presence of long tubular structures in the soil, or "macropores," could significantly 
influence the flux of water through the medium. We have not observed any such pores. 
The material is cohesionless and cohesionless grains cannot sustain voids and therefore 
macropores are not expected. According to the observational and theoretical evidence, 
we discard this possibility (no gophers). Local density vari ations may be possible. 

In a porous medium, the velocity is given by Darcy's law ijPolubarinova-Kochinal 1 962|) 



where k is the seepage constant, p is the pressure, p the density, g Earth's acceleration 
constant, and g the unit vector in the direction of gravity. (The permeability n of the 
medium is directly related to the seepage coefficient through = k/{pg), where p, is 
the viscosity of the fluid.) Since the fluid is incompressible, V ■ v — and it follows 
from ()4.1|l that 



Boundary conditions are shown in Figure [TUl for a stationary situation. The bottom of 
the tank is impermeable, the water column in the back is represented by a linearly varying 
pressure, and atmospheric pressure is applied at the seepage surface. At the water table 
there is atmospheric pressure, p = 0, and the velocity is tangential. The shape of this 
free surface adjusts such that both boundary conditions are satisfied simultaneously. The 
pressure distribution and the shape of the free surface are independent of the parameter 
fc, since neither the Laplace equation (|4.2|l nor the boundary conditions depend on k. 
The shape of the water table is also independent of g, since the boundary conditions 
and the Laplace equation can be written in terms of p/{pg)- (The limit g — *■ poses no 
physical problem, since for 5 = every shape is permissible.) While velocities depend on 
k and g, the water table is purely determined by the geometry of the boundaries. This 
is als o obeyed by the known exact solutions for unconfined flow ( Polubarinova-Kochin^ 



In the above formulation of flow in a porous medium, capillarity and evaporation are 
ignored. Evaporation should be negligible over the time period the channels form. On 
the beach and in the experiment the capillary height exceeds the depth of the chan- 
nels, and hence the capillary forces may prevent the detachment of the water table from 
the surface. However, the capillary effect vanishes for larger scale landscapes. Further 




(4.1) 



V^p = 0. 



(4.2) 



Il962l Chap. VII). 
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Figure 11. Cross section of the beach with water table and velocities obtained by numerical 
solution of the porous media equations. At high tide the sand is soaked with water. As the 
tide recedes, water seeps out along the beach face, starting at the point where the water table 
intersects the surface. The effect of capillarity is neglected. 



consideration of the potential role of capillary forces will be necessary before detailed 
quantitative comparisons with experiments are made. We expect however, that our qual- 
itative conclusions will be independent of these issues, because ground water piracy takes 
place in either case. 

4.2. Numerical Method; Planar Sloped Surface 
IPolubarinova-Kochinal ( 19621) provides a comprehensive collection of analytic solutions to 



problems of groundwater movement, but most geometries need to be solved numerically. 
We solve the Laplace equation numerically with a conventional relaxation method on 
a regular grid. The scheme uses the pressure at the six neighboring grid points and 
iterates until the pressure field has converged. At the free boundary, cells are no longer 
treated as rectangular, but in the three-dimensional finite difference approximation of 
the Laplace operator, the exact location of boundary point(s) is taken into account. 
At the free boundary, atmospheric pressure is imposed and the free surface is moved 
until it is stationary, using a simple first-order time integration. The time-advanced grid 
points representing the free surface no longer lie on a regularly spaced grid and are 
interpolated back onto a rectangular grid at each integration step. In three dimensions, 
when the free surface is two-dimensional, the interpolation fits a plane through three 
nearby points. The three points are chosen such that without variations along the x- 
axis, one exactly recovers a one-dimensional, linear interpolation scheme. Hence, three- 
dimensional simulations without variation along the cc-direction yield exactly the same 
results as two-dimensional simulation s. The numerics are tested by comparison with the 
analytic solution of the vertical dam ( Polubarinova-Kochinal[l962 , eqs. Vn.10.34, 10.35, 



10.41 ) and with analytical estimates for the flux through a wedge ( Polubarinova-Kochinal 
Il962l eq. VII.12.12). 

FiEfurelTTIshows the velocity field in a geometry similar to the beach (without channels). 
As the height of the water table decreases along the horizontal direction, the velocity 
increases as a result of the incompressibility of the fiow. 

Th ere are approximate solutions for a wedge without tailwater ijPolubarinova-Kochinal 
[1962^. If the slope s is small and the wedge is long, then the height ho at which the water 
table intersects the surface (see Figure is approximately h^ k, [Q /k){l/2 + 1/s), 
where Q is the flux per unit width and k the seepage constant. The average horizontal 
velocity at the point where the water table intersects the surface is w = Q/ho, hence 

2ks 



ks. (4.3) 

z, -r s 

The velocity increases with slope, because ho decreases with slope. Since the velocity 
changes little along the vertical (see Figure [TT|l . the average velocity v is a, good approx- 
imation to the velocity at the intersection, vt ~ v. Hence, the velocity with which water 
emerges at the seepage face is higher for higher slopes. Our numerical simulations show 
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that for sfopes of 10% or smaller the presence of tailwater does not change the velocity 
vt significantly. 

4.3. Wentworth Instability 

So far we have only considered seepage without channels. Channels can be represented 
in the simulation by applying the atmospheric boundary condition at a properly shaped 
surface. This shape is the planar seepage surface with channels carved out. The depth of 
water running within channels is neglected. 

The presence of a channel perturbs the flow field. Once a channel is formed, flow lines 
focus toward it, but furthermore it leads to inclination of the underground water table in 
response to an outlet channel. Figure [T^ illustrates the deformation of the underground 
water table. This brings even more water to the channel and promptly acts as an e f ficient 
instability mechanism. This scenario was qualitatively described bv I Went wort hi l)l928() 
in the context of Hawaiian valleys. Hence, we refer to it as "Wentworth Instability". 
Figure shows the velocity field on the water surface. The channel collects flow over 
a certain lateral distance, hence suppressing the growth of other channels nearby, due to 
groundwater competition. 

The area over which the channel draws water can be visualized by following the flow 
lines which reach from the water reservoir to the bottom end of the channel. The vertical 
component of velocity, Vz , is small and approximate flow lines are obtained by integrating 
the horizontal components of the velocity fleld. The geometry of the region of influence 
(the shaded region in Figure 112b ) is again purely determined by the geometry of the 
boundaries, including the dimensions and position of the channel. It is independent of 
density, permeability, and other material properties. The extent of this drainage area (or 
drainage volume) describes the screening effect of the channel and is therefore a relevant 
mathematical problem in this context. 

The overall downhill velocity Vy « vt is approximately constant along the seepage face. 
The width W over which water is collected into a channel of length £ approximately obeys 
W/£ oc Vx/vy. As discussed above, variations of Vx and Vy in space are slow, and hence 
the width W is roughly proportional to channel length £. The constant of proportionality 
depends on channel depth (via Vx) and on the slope (via Uy), unless W extends to a 
distance 0{t) over which Vx changes. Figure El in comparison with Figure [T2b. shows 
this effect in full numerical simulations. The width of area drawn into the channel becomes 
smaller and smaller for shorter and shorter channels. 

Groundwater piracy should relate to channel spacing. Since the lateral claim of chan- 
nels grows proportionally with their length, it is small initially. Small channels have 
only a small influence on the flow fleld. Long channels on the other hand will compete 
with nearby channels. With the growth of the channel, its region of influence on the 
groundwater movement grows. 

4.4. Discussion 

4.4.1. Length-scale for the width of the subsurface drainage area 

In this section we further discuss the width over which channels draw groundwater. 

Figure [HI illustrates, in a two-dimensional auxiliary problem, that the influence of a 
channel of size e is not 0(e), as it would be without a free surface, but 0{L), where L is the 
distance to the boundary. Without free surface, if only the atmospheric pressure boundary 
condition were applied, the influence of the channel would only reach 0(e). (Since the flow 
is driven by height differences, the length-scale of the water table deformation translates 
into the length-scale over which the velocity component, Vx^ changes.) The velocities 
themselves strongly depend on e, but the scale over which they change does not. As a 



Channelization in Permeable Ground 



13 





0.4 



0.3 Z 
0.2 Z 



_ _ _ ^^^^^^^^ 



0.1 - : 



0.8 



0.6 

y (m) 



0.4 



0.2 



Figure 12. a) Deformation of the water table due to an outlet channel, as obtained by numerical 
solution of the porous media equations. The aspect ratio in (a) is distorted to emphasize vertical 
changes. The actual slope of the wedge is 10%. The deformations of the water table increase 
the flux of water into the channel, b) Velocities on the surface of the water table, showing the 
influence of an outlet channel on the flow field. (The resolution of the calculation is higher than 
the number of arrows plotted.) Water within the shaded region flows into the channel. Black 
triangles mark the intersection of the water table with the surface. 
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Figure 13. Compared with Figure [T^ . the region of influence changes with the length of the 
channel. All other parameters are the same as for Figure IT^ 



14 



N. Schorghofer, B. Jensen, A. Kudrolli, and D.H. Rothman 



u 1 

^ 2L ^ ►x 

Figure 14. Deformations of the water table (dashed line) in response to a channel of width 
and depth e. Note the large horizontal extent, 0(L), of the deformation. 

consequence, the lateral influence of the channel can be much larger than its width or 
depth. The deformations of the water table, although small, enhance the lateral influence 
of the channel and lead to quick growth of the channel. 

This length-scale can be argued for as follows. Due to incompressibility Vx{x)h{x) ~ 
constant, where h{x) is the height of the free surface. When the behavior is close to 
hydrostatic, p{x,z) « pg{h{x) — z). The actual depth over which horizontal motion 
takes place does not need to be known, but is assumed roughly constant with x. Using 
equation H4.1|l . the velocity is 

k dp dh{x) 

Vx{x,z) — sa ~k— — . 

pg ox ax 

Therefore, h{x)dh{x)/dx « constant. After integration, the unknown constants can be 
expressed in terms of h{0), the uppermost point, and h{L), the lowest point of the free 
surface. We obtain 

h^ix) = h'^{Q) - [h^{Q) - h^iL)] y. (approximation) 

h{L) cannot be determined within this approximation. After differentiation at x = 0, 

_ hjO) + hjL) hjO) - h{L) _ f h{Q)-h{L) \ 
2h{0) L ~ \ L )■ 

The inclination at the boundary is comparable to the average slope of the water table. 
Hence, the water table is significantly sloped over the entire length, not just near the 
channel. 

Full numerical simulations bear out this effect in three dimensions. Figure [T5k shows 
the pressure in an inclined plane 1cm beneath the surface. (Due to the atmospheric 
pressure boundary condition, there are no pressure variations on the seepage surface 
itself.) The distance over which pressures decay corresponds to the distance over which 
the velocity decays in Figure IT^. For the purpose of illustration. Figure shows 
the (unphysical) result if the lateral deformation of the water table were neglected. The 
lateral influence of the channel is comparable to the channel width, while in Figure [T3k 
it is proportional to the channel length. 

4.4.2. Coarsening 

The theoretical description is incomplete without quantitative relations that describe 
the rate of erosion as a function of the amount of water seeping out the surface and 
pouring down the channels. Assuming that erosion is a monotonically increasing function 
of the amount of water drawn into the channels, it is possible to predict the dynamical 
behavior qualitatively. 

Groundwater piracy should relate to channel spacing. The "wavelength" A ~ c'W. Here 
c' is a proportionality constant c' < 1, since piracy only gradually slows the growth and 
does not instantly set the spacing. According to the previous discussion, W is proportional 
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Figure 15. The pressure distribution 1cm beneath the surface, on a plane parallel to the 
slope. Contour intervals are lOPa. A darker gray tone indicates lower pressure, a) The actual 
pressure distribution, b) Pressure distribution when the lateral deformation of the water table 
is neglected. This case is unphysical, but illustrates that the response of the water table is much 
shorter-reached. In this simulation the overall water table is first calculated without channel. 
Starting with this surface the channel is emplaced and the pressure field calculated. Only in (a) 
are resulting lateral deformations of the water table taken into account. The channel is 1cm deep 
(perpendicular to the x — y plane), 1cm wide (parallel to the x-axis), and 40cm long (parallel 
to the y-axis). 



to £, 

W = c£, (4.4) 

where c is a geometric constant, independent of permeability or erosion threshold. On 
the other hand, c may depend on slope and channel depth. The resulting proposition is 

A«c'c£. (4.5) 

Relation H4.5I) implies that seepage erosion has no constant preferred wavelength, since 
the channel length is time-dependent, i — £{t). This is consistent with the observation 
of coarsening in the experiment. Formula H4.5|l is not applicable to the early stage of 
evolution, which is strongly influenced by the initial incision processes. 

Characteristic channel spacing can simply arise from the presence of a dominant length 
scale in the longitudinal (flow) direction. That is, if the length of the channel is roughly 
predetermined, then the drainage outlets would be spaced at a distance roughly propor- 
tional to the length of the channels. The Wentworth instability predicts that channels 
keep on growing and competing, until the externally given maximum length is reached. 
In a natural setting the maximum length may be determined by the length over which 
erosion occurs. 

The scale-separation in the experiment may not be large enough to distinguish whether 
spacing is proportional to channel width or channel length. However, in both theory and 
experiment competition for ground water controls channel development. The coarsening 
mechanism is qualitatively the same. 

In summary, coarsening of the lateral inter-channel spacing is predicted by theoretical 
considerations. Channels compete for groundwater, which slows and eventually stops the 
growth of some of them. The length-scale characterizing the distance between major 
channels increases with time, until the channels have reached their maximum length. 
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Eventually the lateral spacing should be only confined by the externally constrained 
maximum length. 

5. Conclusions 

The lessons wo have learned arc about the role of runoff and seepage during channel 
formation, the generic nature of rhythmic channelization, and the importance of the 
water table's free surface. 

We observe the formation of drainage outlet channels in permeable ground, on the 
beach as well as in laboratory experiments. In both cases, initial channels are quickly 
incised by surface flow. Even in a very permeable medium, small amounts of surface 
runoff are more efficient in starting channels than is the emerging groundwater. Following 
the incision, groundwater from the vicinity seeps into the channel and flows inside the 
channel, enlarging it, mainly lengthwise, through erosion. 

An initial lateral separation may be related to the initial sheet flow. Coarsening of 
the lateral inter-channel spacing is discernible in the experiment. For narrowly spaced 
channels, one channel often slows and eventually stops its growth, presumably due to 
piracy of groundwater by a nearby channel. A second process taking place is the belated 
formation of new channels far away from existing ones. Rhythmic channel patterns arise 
as water seeps through a mono-disperse permeable material and erodes material as it 
emerges and pours down the channel. 

These observations lead us to a theoretical description in terms of flow of ground water 
in a homogeneous porous medium with a free water table. According to this theory, the 
response width of the water table is much wider than the channel is wide or deep and 
approximately proportional to the channel length. This underscores Wentworth's insight 
that this an efficient (and therefore the dominant) instability mechanism, because it is 
prompt and far-reaching. The groundwater flow pattern, for a given channel geometry, 
is independent of material properties. 
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